Typeset with jpsj2.cls <ver.l.l.l> 



A Toy Model of Flying Snake's Glide 

1 Koji Matsumura* and 1,2 Y-h. TagucW 

1 Department of Physics, Faculty of Science and Technology, Chuo University, 1-13-27 Kasaga, Bunkyo-ku, Tokyo 

112-8551. 

2 Institute for Science and Technology, Chuo University, 1-13-27 Kasaga, Bunkyo-ku, Tokyo 112-8551. 

(Received February 8, 2008) 

We have developed a toy model of flying snake's glide [J.J. Socha, Nature 418 (2002) 
603.] by modifying a model for a falling paper. We have found that asymmetric oscil- 
lation is a key about why snake can glide. Further investigation for snake's glide will 
provide us details about how it can glide without a wing. 
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1. Introduction 

Biological fluid mechanics 1 ' 2 ) is recently an interesting field. It is unclear why dolphins or tuna 
fishes can swim so rapidly and why small insects like flies or bees can get large enough lift force. 
However a flyingsnake 3 ) is stranger than them in the light of fluid mechanics. Although it does 
not have anything like wings, it can glide with constant falling speed. This means it can produce 
at least some lift force which cancels downward gravity. Here we have shown how it can glide 
without wings using the toy model 4 ) which is originally proposed for explaining the behaviour of a 
falling paper. A falling paper exhibits several behaviour ranging from a simple falling to a periodic 
rotation. 5 ^ However, there are no directional gliding motions which the flyingsnake can have. We 
have found asymmetric oscillating motion can induce a gliding motion. Probably this is the reason 
why the flyingsnake must oscillate when it glides. 

When the flying snake would like to glide, 3 ) it jumps up from some higher location. After some 
short transient time period, it starts to glide. While gliding, it shapes on horizontal plane a S letter. 
Furthermore head-to-tail distance oscillates periodically. It seems to get lift force without doing 
anything special other than that. It is the rather difficult task to understand by fluid mechanics this 
mechanism. Fluid around a falling snake is disturbed violently and its spatio-temporal patterns are 
very complicated, probably something like that called turbulence appears. It is difficult not only 
to understand such complicated spatio-temporal patterns but also to simulate them numerically. 
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2. A Falling Paper Model and its Modification for Flying Snake 

Tanabe-Kaneko model 4 ) can reproduce qualitatively the behaviour of a falling paper; chaotic 
or periodic rotations, chaotic or periodic flutterings and a simple falling, although there are some 
criticisms to this model. 6 ** In Tanebe-Kaneko model, a falling paper has only three state variables; 
a horizontal velocity, a vertical velocity, and an angle of inclination of paper. This means, a falling 
paper is modeled as a line segment which goes down in vertical plane. Thus its motion is restricted 
to be within two dimensions. External forces are a gravity, a viscosity force and a lift. Although 
the later two forces are complicated functions of the velocity field of the surrounding fluid, Tanabe 
and Kaneko assume they depend upon only a horizontal velocity, a vertical velocity and an angle 
of inclination of paper. In this study we employ this model in order to describe a flying snake's 
gliding motion. 

In order to model snake's gliding with this model, we have assumed that 

(1) Neither viscosity force nor lift changes due to oscillations. Viscosity force is mainly dependent 
upon level cross-section area on which snake shapes S letter. This area does not change due to 
oscillation. Lift mainly depends upon snake's volume, which is conserved during oscillation. 

(2) Inertia moment is strongly dependent upon oscillation because head-to-tail distance is propor- 
tional to its gyration radius whose squared value is proportional to inertia moment. 

Except for those mentioned above, snake is regarded as a falling paper. 

When snake oscillates, the distance d between head and tail has time dependence as d(t). Other 
than making inertia moment dependent upon time, we do not modify Tanabe-Kaneko model at all. 
As a result, our model for flying snake is 7 ) 

it = — (k± sin 2 9 + k\\ cos 2 9)u + (k± — k\\ ) sin cos 9v =F KpV 2 cos(a + 9) cos a, 

v = (k± — ku ) sin 9 cos 9u — (k± cos 2 9 + fen sin 2 9)v ± irpV 2 cos(a + 9) sin a — g, (1) 

Co = —k±LO — (3irpV 2 /d{t)) cos(a + 9) sin(a + 9), 

where fluid density ratio p is pf£ c /m pi pf is fluid density, m p and l c are mass and horizontal size 
of snake respectively. As mentioned above, l c does not change with oscillation but takes constant 
value, u and v are horizontal velocity and vertical one, ui is the angular velocity on the vertical 
plane. k± (k\\) is the friction coefficients along the perpendicular (parallel) direction to a plane on 
which a snake shapes a S letter. 9 is incline angle of this plane, a = arctan(u/-u) and V = Vu 2 + v 2 . 

3. Results 

We have tried to find d(t) which induces the gliding motion which Tanabe-Kaneko original 
model lacks. First of all, we employ simple harmonic motion as d(t). It cannot induce gliding 
motions at all, although we check any possibilities that original model exhibits any motions. Thus 
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we tried other oscillating motions. After several trials and errors, we found three cases which can 
produce gliding motion. These three d(t)s are 

0.3(1 + sino;ot) 2 + 0.7, 
d(t) = < exp(— 1 — 1.3sina;oi), 
l/[l+0.5sinwo*]- 

In the case of all of these three, d(t) oscillates around 1.0. When d(t) takes the constant values, 
this corresponds to original model of a falling paper. Also, as long as we tried, gliding motion can 
take place only when d(t) oscillates around the value with which a periodic fluttering occurs in the 
original model. In Figure 1, we show periodic fluttering of a falling paper. We also show 9(t) as a 
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Fig. 1. Left: trajectory of periodic fluttering of a falling paper reproduced by Tanabe-Kaneko model (d(t) = 1). 
Right: 0(t) as a function of u(t). Their relation is the same as that between the velocity and the coordinates 
observed when harmonic motion takes place. 

function of u(t). Due to the symmetric motion, their relation is also symmetric and is nearly equal 
to that between the velocity and the coordinates observed when a harmonic motion takes place. 

Now in order to simulate snake case, we have to consider the time dependence of d(t). In Figure 
2(a), we show d(t) as a function of t. This aymmetric oscillation has turned out to be important 
in order to induce gliding motion. This is the reason why simple harmonic motion cannot induce 
gliding motion. Due to this asymmetry, as d(t) becomes larger, velocity of oscillation becomes 
higher, in contrast to the simple harmonic motion. When oscillation is asymmetric as shown 
above, 9{t) can resonate to the oscillation (Figure 2(b)) if ujq is chosen properly. This means that 
the period of fluttering is forced close to that of oscillation. Moreover, as shown in Figure 2(c), this 
resonance violates time reversal symmetry of 9{t) which 9{t) for periodic fluttering has. Because 
of this symmetry breaking, periodic fluttering motion becomes gliding motion although fluttering 
motion does not vanish completely (See Figure 3). For a long period, a falling paper keeps constant 
velocity and constant angle 9(t) as if flying snake keeps them. 3 ) Also in Figure 3, we show 9(t) 
as a function of u(t). Since asymmetric motion of 9(t) results in that of u(t), in average a falling 
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Fig. 2. (a) d(t) — 0.3(1 + sinuiot) 2 + 0.7, loq = 0.9 (b) 0(t) when d(t) = 1. Its time dependence is almost sinusoidal 
because paper flutters periodically, (c) 6(t) when varies as shown in (a). Due to the oscillation, period of 
fluttering becomes close to that of oscillation. 




Fig. 3. Left: A trajectory in vertical plane. Right: 8(t) vs u(t) for the left figure. The symmetry seen in the motion 
without oscillation (Figure 1) has vanishes. 
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paper can move along one direction. This is the reason why oscillation causes gliding motion, i.e. 
a falling along one direction. 




horizontal distance t U 



Fig. 4. Left:d(i) = exp( — 1 — 1.3sinwo£), too = 0.9 Center: A trajectory in vertical plane. Right: 6(t) vs u{t) for the 
center figure. 
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Fig. 5. Left:d(t) = 1/[1 + 0.5sino>n£], wo = 0.9 Center: A trajectory in vertical plane. Right: 6 it) vs u(t) for the 
center figure. 

For the remaining two cases, we also show trajectory, d(t), and 9{t) vs u(t) (Figures 4 and 5). 
In order to induce gliding motion, for these two cases, d(t) must have time dependencies similar 
to that of the first case, too. Although all of them induces the motions along the one direction, 
there are some differences among these three cases. For the first case, 9 — u plot shows very narrow 
band, which means gliding motion is very periodic. On the other hand, for the second case, 9 — u 
plot tends to converge to the periodic motion very slowly. And for the third case, 9 — u plot shows 
relatively broader band, this means simple periodic motion is modulated by slower periodic motion. 

In spite of these minor differences, their global behaviours are very similar. Thus, there should 
exist much more such possibilities of d(t) that cause gliding motion. Of course, d(t) cannot take 
so different shapes from the above cases. As long as we tried, only these shapes can induce gliding 
motion. Probably, a flying snake would have found such very rare possibilities during its evolutions, 
because it has much more time for trial than we spent. 

It is also a problem that we cannot stop fluttering motion completely. However, a flying snake 
can choose any other much more complicated time dependence of d(t) Further detailed investigation 
of snake's motion will tell us which d(t) can induce a simple gliding motion that is not accompanied 
with the fluttering. 
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4. Sensitivities to the detailed parameter values or the forms of the function d(t) 



Finally, we would like to comment on sensitivities to the detailed parameter values or the forms 
of the function d(t). Once we fix all parameter values other than those included in d(t), we can 
change the functional form of d(t) very little. For example, if we slightly change the value 0.3 in the 
form of d(t) = 0.3(1 + sin wot) 2 + 0.7, the directional gliding motion disappears quickly. Actually 
speaking, three d(t)s mentioned above are tuned so that they have almost equal time dependency. 
Thus, we believe that this specific form is very important. 

Also, the condition for the directional gliding motion is asymmetry. As mentioned above, 
asymmetry in d(t) results in that in velocity which causes directional motions. For example, if we 
employ a normal sinusoidal motion 

d(t) = 0.6 sin u t + 1.3, u = 0.9, 



there is not any asymmetry in velocity, thus directional gliding motions cannot appear (Fig. 6). 
In spite of the fact that d(t) is tuned so that its functional form is close to those of three d(t)s 




Fig. 6. Left:d(i) = 0.6 sin ui t + 1.3, wo = 0.9 Center: A trajectory in vertical plane. Right: 9(t) vs u(t) for the center 
figure. 

excluding the asymmetry, directional motion cannot occur. Thus it is clear that asymmetry is the 
most important factor. 

It is also easily understood that the direction of motion is determined by the initial condition. 
Asymmetry in d(t) cannot decide in which direction symmetry is broken. This depends upon how 
the snake starts to fly. 
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